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ABSTRACT 

We report the results of a survey for fluorescent Lya emission carried out in the field surrounding 
the z = 3.1 quasar QSO0420-388 using the Focal Reducer/Low Disperion Spectrograph 2 (F0RS2) 
instrument on the Very Large Telescope (VLT). We first review the properties expected for fluorescent 
Lya emitters, compared with those of other non-fluorescent Lya emitters. Our observational search 
detected 13 Lya sources sparsely sampling a volume of ~ 14000 comoving Mpc'^ around the quasar. 
The properties of these in terms of (1) the line equivalent width, (2) the line profile and (3) the value 
of the surface brightness related to the distance from the quasar, all suggest that several of these may 
be plausibly fluorescent. Moreover, their number is in good agreement with the expectation from 
theoretical models. One of the best candidates for fluorescence is sufficiently far behind QSO0420-388 
that it would imply that the quasar has been active for (at least) ^--^60 Myr. Further studies on such 
objects will give information about proto-galactic clouds and on the radiative history (and beaming) 
of the high-redshift quasars. 

Subject headings: cosmology: observations — intergalactic medium — quasars: individual: QSO0420- 
388 — line: identification — galaxies: high-redshift 



1. INTRODUCTION 

The analysis of absorption systems in the spectra of 
high-redshift quasars has represented for several decades 
the only practical way to obtain information about the 
properties of the intergalactic medium (IGM). The Lya 
forest gives an unique insight into the one-dimensional 
distribution of low density hydrogen along the quasars' 
line of sight. Higher column-density features, i.e. the 
Lyman-limit systems (LLSs) and the damped Lya sys- 
tems (DLAs), correspond to denser concentrations of 
atomic hydrogen and, for this reason, can be associated 
with galaxy formation. Unfortunately the information 
from quasar absorption spectra is almost always one- 
dimensional. 

An interesting alternative to absorption studies is to 
try to detect the IGM in emission rather than in absorp- 
tion. The absorption of ionizing photons should be as- 
sociated with the emission of fluorescent Lya photons, 
as originally proposed by Hogan & Weymann (1987). 
Detection of the fluorescent emission would provide a 
three-dimensional picture of the neutral IGM and proto- 
galactic clouds at high redshift. 

In an earlier paper (Cantalupo et al. 2005, hereafter 
Paper I), we constructed models of fluorescent emission 
in realistic hydrodynamical simulations. These simula- 
tions were normalized so as to produce the correct statis- 
tics for Lyman limit absorption systems. For clouds that 
are optically thick to the ionizing radiation (i.e. Ly- 
man limit systems and above), the surface brightness 
(SB) of the fluorescent emission is set by the strength 
of the ionizing background. Unfortunately, the inten- 
sity of the UV background at z ~ 3 (e.g. Haardt & 
Madau 1996) corresponds to a maximum Lya SB of 
~ 3 X 10~^°erg s~-'^cm~^arcsec"^ . Detecting such a low 
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SB is very challenging with present-day instrumentation. 
Blind searches in the held have only produced a number 
of null results (Lowenthal et al. 1990; Martmcz-Gonzalez 
et al. 1995; Bunker, Marleau & Graham 1998). 

The SB of fluorescent emission will be higher if the lo- 
cal ionizing background is increased, e.g. for gas clouds 
lying close to a bright quasar. In Paper I, we examined 
this effect and showed that it could boost the fluores- 
cent emission to detectable levels over a large volume, 
i.e. of order 30,000 comoving Mpc'^ for a sufficiently 
bright quasar. The penalty for this boost is that the ex- 
tra photo-ionization increases the total hydrogen column 
density that is required to reach a given neutral column 
density, i.e. to produce the maximum fiuorescent surface 
brightness. This is the analogue of the proximity effect 
seen in quasar absorption spectra. It means that maxi- 
mally fluorescent sources near to a quasar will be systems 
that, in the absence of the quasar, would be DLAs rather 
than LLS. DLAs have hydrogen column densities compa- 
rable to a sightline through the disk of present day spiral 
galaxies and a total gas mass similar to the mass in stars 
and gas within galaxies at the present epoch (Wolfe et 
al. 1995). For these reasons, DLAs are widely believed 
to be the gas reservoirs from which the present galaxies 
formed. However, it is still unclear if DLAs are galaxies 
already in place (e.g. Prochaska & Wolfe 1998) or are 
still in the form of protogalactic clouds (e.g. Haehnelt, 
Steinmetz & Ranch 1998). They typically contain very 
little of the molecular hydrogen usually associated with 
star formation (e.g. Ledoux et al. 2003, but see also 
Zwaan & Prochaska 2006). Moreover, metallicity mea- 
surements show a large range from 1/10 solar to a value 
close to that measured for the Lya forest (cf. Prochaska 
et al. 2003 and Simcoe et al. 2004). It is likely that 
DLAs include systems of different origins. 

In addition to studying the intergalactic medium, the 
detection and measurement of fluorescent Lya emission 
from clouds at different distances and orientations from 
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the quasar, would in principle allow the study of the his- 
tory and angular distribution of the ultraviolet emission 
from quasars. Both questions are of considerable topical 
interest in the context of unified models of active galactic 
nuclei (see e.g. Antonucci 1993). 

In this paper we describe a blind survey for fluo- 
rescent Lya emission carried out around the z ^ 3.1 
quasar QSO0420-388, one of the brightest quasars at high 
redshift for which the Lya line lies within an existing 
narrow-band filter of the VLT-F0RS2 instrument. 

Previous searches for fluorescence around quasars have 
had either null results at z ^ 2.2 (Francis & Bland- 
Hawthorn 2004) or only one (or marginally two) pos- 
sible detection at 2 ^ 4.3 (Francis & McDonnell 2006). 
In the first case, the null result is perfectly consistent 
with our model of fluorescent Lya emission presented in 
Paper I. In the second case, the very close association 
of the system with the quasar itself makes the identifi- 
cation with fluorescence difficult (see section 2.5 for de- 
tails). Other possible fluorescent sources have been de- 
tected serendipitously in connection with absorbing sys- 
tems in quasar pairs (Fynbo et al. 1999, Adelberger et 
al. 2006). However, the metallicity of the systems and 
the presence of a substantial continuum emission prob- 
ably indicate the presence of star formation within the 
clouds. Lya emission at high redshift can of course be 
produced by internally ionized clouds, either from an as- 
sociated young stellar population or from processes in- 
volving the gas itself (e.g. cooling). Several observations 
during the last decade have shown that these emitters are 
quite common in the high-redshift universe (for a review 
see Taniguchi et al. 2003) and enhanced by clustering ef- 
fects around bright sources like AGNs and, in particular, 
radio-galaxies (e.g. Vcnemans et al. 2005 and references 
therein). In searching for fluorescence, we must thus take 
into account the presence of this population. 

The layout of the paper is as follows. In §2 we com- 
pare the characteristics of fluorescent Lya emission with 
those of other sources of Lya emission at high redshift, 
in order to determine which properties can best be used 
to recognize fluorescent emission. Although fluorescent 
sources should have characteristic equivalent widths (ex- 
tremely high), SB and spectral profiles, none of these 
provide a unique diagnostic. We conclude this section 
by reviewing recent claims in the literature for detected 
fluorescence. In §3 we describe our own new observa- 
tions and the data reduction. We then isolate a set of 
Lya emitting candidates around the target quasar. In 
§4, we discuss in detail the properties of the individual 
sources in the context of the expected fluorescent char- 
acteristics, and compare their number density with both 
our own models for fluorescence (Paper I) and those of 
other non-fluorescent Lya sources. We conclude that it is 
likely, but not certain, that some of the detected sources 
that we have detected are indeed fluorescent. In §5 we 
summarize the paper. 

2. RECOGNITION OF FLUORESCENT EMISSION 

Besides fluorcsccint emission, Lya can be produced by 
several kinds of internal sources, e.g. photo-ionization 
from a young stellar population or active nucleus (LAE, 
Lya emitter [galaxies]), or from cooling of shock-heated 
gas. Many observations during the last decade have 
shown that Lya emitters are quite common in the high- 



redshift universe (e.g. Cowie & Hu 1998: Pascarelle, 
Windhorst & KeU 1998; Fynbo et al. 2003) and enhanced 
by clustering effects around objects like AGNs and, in 
particular, radio-galaxies (e.g. Le Fevre et al.l996; Pas- 
carelle et al. 1996; Venemans et al. 2005) which likely 
reside in high density regions of the Universe. 

In this section, we review the characteristic properties 
that in principle could identify a source as fluorescent: i) 
the equivalent width (EW), ii) the emission profile, iii) 
the surface brightness (SB), and iv) the number density 
of emitters. We will take advantage of the model re- 
cently presented in Paper I, in which a combination of a 
high-resolution hydrodynamic simulation and two radia- 
tive transfer schemes was used to predict the properties, 
spectra and number density of fiuorescent Lya emitters 
at 2: 3. 

2.1. Equivalent Widths 

In the fluorescent case, the Lya emission origi- 
nates from pure gaseous clouds externally photo-ionized. 
Therefore, any continuum emission comes only from re- 
combination processes: mainly two-photon continuum 
and free-bound recombinations. However, given the low 
value of the emission coefficient corresponding to these 
processes, the continuum level around the Lya wave- 
length is very faint, resulting in an extremely high EW. 
The case of Lya emission originating from cooling shock- 
heated gas will be similar. 

If the Lya emission results from photo-ionization by 
embedded young stars, we expect a significant contribu- 
tion to the UV continuum emission from the stellar pop- 
ulation. In this case, the Lya EW depends on the age, 
initial mass function (IMF), metallicity and dust con- 
tent of the star-formation burst. Synthesis population 
models (see e.g. Schaerer 2003 and references therein) 
show that high EWs can be produced for very young 
bursts in metal-poor and dust-free galaxies. Upper limits 
for very metal-poor models {0 < Z < 10"^) range from 
EW~ 8OOA (for a very young population with an age of 
2 Myr) to EW-- 200A (for a 10 Myr old star-formation 
burst). Synthesis models of star- forming galaxies with 
moderate ages (> 10 Myr) and with metallicity (and 
IMF) closer to the normally observed values have EW in 
the range 50 - 200 A (Chariot & Fall 1993). Surveys of 
Lya emitters at z ^ 3 usually show a distribution of EW 
in the range from 15 — lOOA, with a few objects extend- 
ing to higher values (up to 200-300A). High-EW objects 
can also be associated with AGN activity. 

2.2. Emission line profiles 

The properties of the Lya line have been widely studied 
analytically under certain approximations (e.g. Neufeld 
1990 and references therein). In particular, for an ex- 
tremely opaque and plane-parallel medium, it is pos- 
sible to demonstrate that the emerging profile consists 
of two symmetric peaks separated by a velocity shift of 
A?; ~ 2(lnrLya)^/^\/2cTth (where Thya and ath are, re- 
spectively, the central-line optical depth of Lya photons 
and the thermal velocity dispersion of the medium). 

Since fluorescent photons originate from externally 
ionized clouds, their typical TLya should correspond to 
TLya ~ 2 X 10^ (assuming a temperature of T = 2 x 10^ 
K). Therefore, in the above approximations, fluorescent 
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emission should be characterized by two symmetric peaks 
separated by Av ^ 8.9(Tth- 

However, detailed simulations of fluorescent emission 
in more realistic situations show that the inclusion of 
velocity fields inside and around the clouds influences the 
emerging line profile. In most cases, the symmetry of the 
double-humped profile is lost and one of the two peaks is 
severely suppressed (see Fig. 7 in Paper I; notice that the 
x-scale of that figure is erroneously multiplied by a factor 
of 4). Furthermore, if the gas consists of several small 
clumps, a turbulent velocity term should be added to the 
thermal velocity dispersion, increasing the separation of 
the peaks. On the other hand, if the turbulent velocity 
of the emitting region is large compared with that of the 
intervening layers, the symmetry of the line profile may 
be restored. 

While the double peaked nature of fluorescent emis- 
sion can be modified, the single Lya lines generated by 
internal photo-ionization at low optical depth can be par- 
tially absorbed by infalling or expanding gas. As a result, 
the emergent proflle can again consist of multiple com- 
ponents. 

We conclude that these practicalities mean that the 
analysis of the emission line profile alone cannot provide 
a watertight signature of fluorescence. 

2.3. Surface Brightness 

With a simplified slab geometry (e.g. Hogan & Wey- 
mann 1987, Gould & Weinberg 1996) of a semi-infinite 
and static medium, fluorescent (self-shielded) clouds 
should act like "mirrors" : illuminated by an external 
source they return the incident ionizing photons in the 
form of Lya line emission with an efficiency of ~ 65%. 
In this case, the Lya SB is uniform and its value is easily 
predicted once the impinging ionizing flux is known, and 
vice- versa. 

In the case of a locally enhanced ionizing flux, e.g. 
from a nearby quasar, the fluorescent SB would be ex- 
pected in this simplified geometry to be higher. We 

introduce a "boost factor" b, which is defined as the 
increase in the predicted SB relative to that induced 
by the UV backgroimd at z ^ 3 (SBbg — 3.67 x 
10-20erg s-icm-2arcsec-2), i.e. SB= (1 -h 6) SBbg. The 
value of b corresponds to: 

^^^g^ i^LL 0.7 / r 

1030 erg s-iHz"^ a VlMpc/ ^ ' 

at a physical distance r from a quasar with monochro- 
matic luminosity Li, = I/LL(i^/t'LL)~" (see figs. 6, 9 and 
section 3.3 in Paper I for further details). 

When velocity fields and a more realistic gas distribu- 
tion are taken into account, it is found (Paper I) that 
the SB may be reduced. In particular, for anisotropic 
illumination (i.e. fiuorescence induced by a quasar), the 
expected SB in the direction of the incoming flux can 
be much lower than expected in the slab approxima- 
tion. The importance of this effect, which we called 
the "geometric effect" . dc;pends on the relative size of 
the shielding- layer with respect to the cloud radius, and 
therefore on the gas density profile and the intensity 
of the impinging radiation. In Paper I, we found that 
the fluorescent SB from self-shielded clouds is still pre- 
dictable in good approximation from the empirical rela- 



tion: 

SB = (0.74 + 0.50 6°-^9)SBbg (2) 

where b represents the ideal "boost factor" described 
above. 

The SB predicted in this way represents a fundamen- 
tal upper-limit constraint: emitters with higher SB must 
receive a ionizing contribution from other sources and, 
therefore, are with high probability internally ionized 
clouds. However, several uncertainties remain to be 
taken into account. In particular, the exact distance 
from the quasar may be uncertain due to peculiar mo- 
tions (and uncertainty in the precise systemic redshift of 
the quasar) . Anisotropy or time- variability in the quasar 
emission will also modify the SB. In fact this latter ef- 
fect offers the possibility, if other indicators indicate a 
fluorescent origin of the emission, that their SBes can 
be used to study the angular distribution and temporal 
history of the quasar UV emission. 

2.4. The number density of emitters 

The covering factor of fluorescent sources on the sky 
within a volume is directly related to the number density 
of optically thick absorbers seen in quasar spectra. This 
fact was used in Paper I to normalise our simulations. 
Two complications enter in the present case. First, the 
regions around quasars are likely to have higher than av- 
erage density. Second, the enhanced ionizing radiation 
will produce a higher photo-ionization in a given cloud, 
reducing the column density of neutral Hydrogen. The 
high surface brightness fluorescence will still come from 
systems which would be seen as LLS if there was a back- 
ground quasar, but these are systems that would have 
been seen as DLA if the extra ionizing radiation had 
not been present. The fluorescence near to a quasar is 
thus closely linked to the "proximity effect" (e.g. Bajtlik, 
Duncan & Ostriker 1988). Since statistics of LLS near 
QSOs are not constrained, we must rely on theoretical 
models to determine their covering factor. 

The procedure to calculate the expected number of flu- 
orescent emitters around a quasar, outlined in section 3.4 
in Paper I, consists of three steps: i) convert the SB sen- 
sitivity limit of the survey (SBum) into a minimum boost 
factor inverting Eq.(2), ii) use b = 6min in Figure 11 
of Paper I to determine the physical number density of 
fluorescent sources with SB>SBiii„; iii) given the values 
of Lll, for the selected quasar, invert Eq. 1 using b = 
bmin to flnd the maximum distance (rmax), and thus the 
volume, where the fluorescent SB>SBiim. As an exam- 
ple, a sensitivity of SBum ~ 10~^^erg s~^cm~^arcsec~^ 
corresponds to a fluorescent (physical) number density of 
dN/dV 0.3 ± 0.1 Mpc~3 (for sources with area A ~ 4 
arcsec^), neglecting any density enhancement around the 
quasar. 

This may be compared with the number density of 
other Lya sources around 2: ~ 3 quasars. Recent surveys 

indicate that the number of Lya emitters around radio- 
galaxies at that redshift can be enhanced by a factor ^ 3 
with respect to the field. In particular, Venemans et al. 
(2005; hereafter V05) found an emitter number density 
of 0.26 ± 0.04 Mpc"^ around the ^ ~ 3.13 radio-galaxy 
MRC0316-257, in a survey with sensitivity of the order 
of SBiim ~ lO^^^erg s^^cm^^arcsec"^ . In the view of 
the AGN unified model (see e.g. Antonucci 1993), the 
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TABLE 1 

Summary of the observations and setups. 



Setup name Date Mode Filter Mask rotation Combined seeing Exposure time 



51 Nov. 29 MXU OIII/3000+51 0° 1".0 25 200 s 

52 Nov. 30 and Dec. 1 MXU OIII/3000+51 90° 1".0 34 200 s 
IMA Dec. 1 Imaging Bcsscl V - 1" .0 7 200 s 

53 Dec. 2 MXU OIII+50 0° 0".7 23 400 s 



environment of radio-loud quasars should be similar to 

that of radio-galaxies. Moreover, we do not expect a 
significant contribution from fluorescent sources to the 
number density found by V05, given the shape of their 
survey and the direction of the (probable) UV emission- 
cone from MRC0316-257. 

2.5. Fluorescent candidates in the literature 

In this section, we briefly describe previous surveys for 
fluorescent emission and review the results of these in the 
light of our discussion above. Given the difficulties in the 
detection of fluorescence induced by the UV-background 
alone, previous surveys have centered around quasars, as 
in our case. 

Francis & Bland-Hawtorn (2004) reported an attempt 
to find fluorescence around a quasar at redshift ~ 2.17 
basing their expectations on the simplifled analytical 
models (Gould & Weinberg 1996). In particular, they 
were expecting to flnd at least 6 fluorescent emitters 
within their sampled region of 1 Mpc around the quasar, 
but they saw none. However, this result is perfectly con- 
sistent with the expectation of our more realistic model 
of fluorescent emission (see section 3.5 in Paper I for de- 
tails). 

Adelberger et al. (2006) reported the serendipitous dis- 
covery of a possible fluorescent emitter associated with 
a DLA at redshift ^ 2.8. In particular they fotmd the 
emitter in the spectrum of a Lyman-break object located 
49" from a much brighter quasar. The object consists of 
two peaks with separation AV ^ 500 km and its SB 
seems to be compatible with the fluorescence induced by 
the quasar, even though the required boost factor is ex- 
tremely high (6 ~ 2000, in the ideal case of a slab). How- 
ever, the object is clearly detected in their G and R im- 
ages (with an observed magnitude Gab = 26.8 ± 0.2) re- 
sulting in a rest-frame equivalent width EWq = 72 ±20 A. 
Therefore, this candidate is more plausibly an internally 
ionized source. A similar detection, associated with a 
lower redshift DLA system {z ~ 1.9) in the spectrum of 
a close quasar-pair, was reported by Fynbo et al. (1999). 
In this case a continuum counterpart is not clearly de- 
tected, although the subtraction of the PSF of the bright 
quasars makes this difiicult. However, the extreme prox- 
imity of the quasar responsible of the possible fluores- 
cent emission (~ 20kpc), indicates a plausible physical 
connection between the two objects. 

Finally, Francis & McDonnell (2006) claimed the de- 
tection of a fluorescent object located just 0.8" from a 
quasar at z = 4.28. Their object is spectrally and spa- 
tially unresolved and the flux is poorly constrained be- 
cause of the proximity of the quasar. There is a lower 
limit for the equivalent width of EWq > 19A. Given the 



vicinity of the quasar (~ 5 — 50 physical kpc) and the 

uncertainties in the flux and size measurement, compar- 
ison with the theoretical SB is hard, although they note 
that the predicted luminosity is ~ 400 times greater than 
observed, if the object is fluorescent. The sampled area 
around the quasar is very small (5" x 7" in projection) 
because of the Integral Field Spectrometer used, making 
any statistical comparison difficult. 

3. OBSERVATIONS AND DATA ANALYSIS 

In this section we describe our own survey for fluores- 
cent emission around QSO0420-388, including the data 
reduction and the selection of the Lya candidates. 

3.1. VLT spectroscopy and imaging 

Observations were taken during four visitor-mode 
nights on the VLT 8.2 m telescope Antu on November 
28 - December 1, 2005. The F0RS2 spectrograph was 
used in multi-object spectroscopy mode (MXU) at stan- 
dard resolution (0.25' x 0.25' arcsec^ pixel"-'^) to build a 
"multi-slit plus filter" (Crampton & Lilly 1999) configu- 
ration centered on the z = 3.1 QSO [HB89]0420-388. 

For all the observations, we used the 1400V grism giv- 
ing a dispersion of 0.61 A pixel"^, in combination with 
the narrow band filters OIII+50 and OIII/3000+51. The 
resolution (R=2100, corresponding to ^ 150kms~^) is 
large enough to resolve the doublet [OII]AA3726, 3729, 
which helps to distinguish a Lya emission from low red- 
shift contaminants. The use of the narrow band filters 
allowed us to put 14 slits (each 6.8' long and 2" wide) 
parallel to each other and separated by 25" without a 
significant overlapping of the spectra on the CCD. The 
total area covered by the slits in the mask corresponds to 
3.17 arcmin^. The same mask was used for all the obser- 
vations, but with two diflferent filters and two different 
spatial orientations: i) slits parallel to North-South di- 
rection, ii) slits parallel to East- West direction. 

The use of two different narrow band fllters and two 
mask orientations was chosen in order to increase the 
sampled volume in redshift and projected space. For our 
setup, the central wavelength of the filters OIII-I-50 and 
OIII/3000-h51 is, respectively 500lA and 5045A with a 
FWHM of 57A and 59A, corresponding to a total redshift 
range z = 3.089 - 3.173 for Lya. Filter OIII/3000-F51 
was used for both the mask orientations (0° and 90° ro- 
tations) and filter OIII+50 for just one orientation (0° 
rotation). 

For each configuration, a series of several dithered ex- 
posures of 1800s each, stepped along the slit by 5 arc- 
sec, were obtained for a combined exposure time of 7 
hours (filter 0III-I-51 and 0° rotation), 9.5 hours (filter 
OIII-h50 and 0° rotation) and 6.5 hours (filter 0III-I-51 
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and 90° rotation). On the third night we also obtained 
a V-band image of the same field for a total exposure 
time of 2 hours. It consisted of 120 dithered exposures 
of 60 sees each, and it was obtained under photometric 
sky conditions. An overview of the observations and of 
the different configurations is summarized in Table 1. 

3.2. Data reduction 

The spectroscopic and photometric data were reduced 
using standard packages (IRAF - Image Reduction and 
Analysis) and additional IDL routines written by the au- 
thors. 

The data reduction consisted of several steps includ- 
ing bias subtraction, spectral distortion correction, fiat- 
fielding (using "sky-flats" obtained from the scientifical 
images themselves) and sky subtraction. The individual 
sky-subtracted images were then combined to obtain the 
three scientific images corresponding to the three differ- 
ent configurations (see Table 1). The combination of the 
images was performed with an averaged sigma clipping 
algorithm in order to minimize the effect of bad pixels 
and cosmic rays. A pixcl-by-pixel noise map associated 
with each science image was generated from the statistics 
of the pixels in the multiple image combinations. 

Spectro-photometric calibration was based on observa- 
tions of different spectro-photometric standard stars for 
each night of observation. The magnitude zero-points 
derived from several standard stars were consistent with 
each other within 0.02 magnitudes. A small Galactic ex- 
tinction correction of Ax{V) = 0.079 mag (Schlegel et al. 
1998) was applied. 

The result of the data reduction was three final "spec- 
troscopic" multi-slit images (SI, S2 and S3; correspond- 
ing to the different configurations of mask orientation 
and filters summarized in Table 1) plus the V-band im- 
age (IMA) of the same field. The ensuing steps included: 
i) the detection of line emission in SI, S2 and S3; ii) 
the identification of possible continuum counterparts of 
the line emission in the V-band image; iii) the selection 
of probable Lya emission with a criterion based on the 
rest-frame equivalent width of the line. 

3.3. Line emission detection 

The detection of line emission candidates was per- 
formed with the help of the software package SExtractor 
(Bertin & Arnouts 1996). In order to make the detec- 
tion process more efficient, we first removed from each 
image the regions corresponding to the brightest contin- 
uum spectra (typically stars and local bright galaxies). 
Such regions were automatically identified for each slit by 
summing the flux over the spectral band and applying a 
sigma-clipping algorithm to the resulting 1-dimensional 
array. We removed all the regions where the integrated 
signal was greater than 3 times the standard deviation 
for that slit. The excluded regions were, respectively, 
~ 4% and ~ 9% of the total image for the 0° and 90° 
rotation configurations. 

The resulting images and corresponding noise-maps 
(obtained in the reduction process described above) were 
then processed by SExtractor. In order to allow the most 
general search, we varied the three most relevant param- 
eters, i.e. the smoothing filter, the minimum detection 
threshold and the minimum detection area (i.e. the num- 
ber of connected pixels). We used always standard filters 



(Gaussian, top hat) plus our own filters with a shape de- 
signed to increase the S/N of the double peaked emis- 
sion. The minimum detection threshold and the number 
of connected pixel was varied, respectively, from 0.3 to 
1.5, and from 3 to 30 pixels. Each combination of the 
three detection parameters produced a separate catalog 
of line emission candidates, for a total of 1344 catalogs for 
each dataset. We refer to these as "positive" catalogues. 

We then applied exactly the same procedure on the 
original images after they had been multiplied by —1 to 
produce a negative image. Any emission objects detected 
in these "negative" catalogues must of course be spuri- 
ous. We therefore kept only those "positive" catalogs 
for which the corresponding "negative" catalogue (gen- 
erated with the same SExtractor parameters) contained 
no detected sources. Notice that this procedure assumes 
Gaussian (rather than Poisson) noise. This is a good 
approximation since the noise is sky-dominated and the 
sky level is high. The final step consisted of merging to- 
gether the lists of objects in these surviving catalogues, 
removing the duplications. 

The output we had from SExtractor, at this stage was 
the central position and approximate size of the can- 
didate line-emission for each of the three spectroscopic 
images. These were used used as a first guess of the 
photometric aperture in the calculation of the line and 
continuum flux (as explained in the next sections). 

3.4. Identification of possible continuum counterparts 

For each candidate detected in the spectroscopic im- 
ages, we identified the corresponding region in the V- 
band image and calculated the continuum fiux fy. The 
size of the photometric aperture was fixed at the slit 
width (2") perpendicular to the slit, and the size of the 
detected object in the spectrum, in the direction along 
the slit. The location of these apertures was obtained 
from accurate relative astrometry of the brightest con- 
tinuum objects which were detected both in the spectra 
and the V image. The errors on fy were estimated from 
integrating on random positions within the image. The 
resulting la limiting magnitude for a typical aperture of 
2" X 2" was about 27.7 mag. 

3.5. Selection of Lya candidates 

In order to select candidate Lya emitters we used a 
method based on the line rest-frame equivalent width : 



EWo = 



line 



Cline(l + Z) 



(3) 



where Fune is the flux of the emission line, Ciine is the 
continuum, if present, at the same wavelength of the line 
and z is the redshift of the emitter. Since the value of 
Ciino is too faint to be reliably measured from the spec- 
tra, we base it on the V-band flux density fy assuming 
a power law continuum with slope /3. Including the con- 
tribution of the line-emission, the fiux density measured 
in the V-band is, to first approximation: 



/v 



-^line • ev (Aline) 

/ev(A)dA 



line 



A, 



eff.V 



A 



line 



(4) 



where ey, and Aeff,v are, respectively, the efficiency and 
the central wavelength (5561. 9A) of the V-band filter. 
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TABLE 2 

Position and redshift of the Lva candidates and the QSO. 



Object 



Position 





'^J2000 


'^J2000 




1 


04:22:18.53 


-38:45:15.6 


3.1197 


2 


04:22:21.17 


-38:46:29.6 


3.1047 


3 


04:22:27.56 


-38:43:24.8 


3.1087 


4 


04:22:27.55 


-38:42:14.5 


3.1599 


5 


04:22:23.36 


-38:45:32.1 


3.1538 


6 


04:22:12.68 


-38:46:39.9 


3.1634 


7 


04:22:19.31 


-38:43:35.8 


3.1298 


8 


04:22:04.06 


-38:45:52.1 


3.1187 


9 


04:22:01.97 


-38:45:59.2 


3.0896 


10 


04:22:23.35 


-38:44:42.0 


3.1263 


11 


04:21:59.87 


-38:43:21.0 


3.0936 


12 


04:22:16.87 


-38:42:27.6 


3.1127 


13 


04:22:21.08 


-38:43:35.8 


3.1298 


QSO 


04:22:14.76 


-38:44:50.7 


3.110'' 
3.123"= 



^ typical error in the redshift estimation is 0.0015 
^ estimated by Osmer et al. 1994 
"= estimated by Lanzetta et al. 1993 



From Eq. (4) we derive: 

.^line 



EWo = 



/v — /line,V \ Al 



Aeff,V 



(l + z)-' (5) 



where 



/line,V 



-Pline • ey (Aline) 

/ev(A)dA 



(6) 



represents the equivalent flux density in the V image 
given by the line emission alone. 

The line flux (i^une) is calculated within an aperture on 
the spectrum (as explained in §3.3). Instead, the value of 
ev (Aline) is interpolated from the measured transmission 
curve of the V-band filter. We have not tried to take 
into account the absorption of the continuum bluewards 
of Aline, because Lya always falls in the blue wing of 
the V-band. We cannot estimate the value of the slope 
(/3) for each source, and so a flat spectrum (/3 = —2) 
was assumed for all the objects. The actual value of /3 
has little effect on EWq - i.e. changing to (3 = —1 (or 
f3 = -3) changes the EW by only ~ 10%. 

In order to compare our results with the other surveys 
for Lja emitters, we selected objects with EWq > IsA 
as candidate Lya line (sec e.g. Vcnemans et al. 2005 and 
reference therein). For candidates with fy < /v(la-) we 
obtain a lower limit to EWq by assuming fy = fy{la). 

3.6. Results 

Each candidate lying above the equivalent width 
threshold (assuming the redshift given from the Lya line) 
was visually inspected in order to remove spurious ob- 
jects like leftover cosmic rays. This resulted in a list of 14 
candidates, of which 12 are compatible to be Lya lines 
while the remaining 2 show a double peaked emission 
that can also be compatible with the [Oil] doublet. One 
of these two candidates has a 2a detection in the V-band 
image and was rejected as a probable Oil emitter. The 
other one, labeled #13, does not show any continuum 
counterpart above la (corresponding to an equivalent 



width greater than 70A for Oil) and, for this reason is 
included in the Lya sample. One other Oil emitter, for 
which the identification is considered secure (see caption 
of Fig. 2 for more details), was found by visual inspection 
of the regions containing a bright continuum that were 
excluded from the analysis at the beginning (see section 
3.3). 

The fraction of contaminants in our sample is thus 
1/14 ~ 7%, similar to the fraction (~ 6.5%) of low red- 
shift interlopers found in other surveys of Lya emitters 
at z = 3.1 (e.g. Steidel et al. 2000; Venemans et al. 
2005). In the following section we will discuss in detail 
the main characteristics of the 13 Lya candidates and, 
in §5, their possible fluorescent origin. 

4. PROPERTIES OF LYa CANDIDATES 

The spatial distribution and the main spectro- 
photometric properties of the 13 Lya candidates are sum- 
marized in Tables 2 and 3. The 2-dimensional spectra in 
S/N units together with the corresponding slit positions 
in the V-band image are shown in Fig. 1, while the 
flux calibrated 1-dimensional spectra are shown in Fig. 
3. The photometric apertures used to calculate the line 
and continuum fluxes are represented by the white boxes 
in Fig. 1. 

Only two objects (#3 and #6) have an integrated con- 
tinuum flux in the V-band image that is well above la 
for the corresponding photometric aperture. However, in 
only one of these (#6) there is a source clearly detected 
at the expected position. In other cases (i.e. #1, #3 and 
#5), the presence of close objects in the T^-band image 
may have affected the; photometry. For the objects with 
fy < la, we estimate a lower limit for EWq, assuming 
fy = la. Thus, only one of our candidates is clearly 
ruled out as a fluorescent emitter based on its EWq. 

The spectral shape of the emission lines appears to be 
sharp and, in most of the cases, consists of only one clear 
peak. One candidate (#4) clearly shows the presence of 
two peaks separated by ~ SA. Both peaks were indepen- 
dently detected by SExtractor at S/N > 5. Candidates 
#3 and #6 also show a lower S/N (> 3) second peak at 
a slightly different spatial position. Candidate #13 was 
detected as a single emission line but subsequent study 
of the unsmoothed spectra revealed the presence of two 
emission lines separated by 3 — 4A. In other words, about 
a third of the sample show some evidence for double line 
structure. 

The estimation of the SB for Lya candidates is compli- 
cated by their small sizes, which typically range between 
the seeing disk and 3 arcsec. For this reason, we have 
simply used the flux in the peak pixel, divided by the 
product of the pixel and the slit size (i.e. 0".25 x 2") 
as a rough estimation of the SB, that will, for the unre- 
solved objects, be a lower limit. Object #4 shows a clear 
flat-topped profile over about 1".25, and for this object 
we estimate the SB directly. 

4.1. Note on candidate #1 

Candidate #1 was detected in a part of the spectro- 
gram in which the spectra of two slits overlap. It was 
assigned to the blue edge of the appropriate slit because 
it lies closer to its central wavelength. We cannot com- 
pletely rule out the possibility that it belongs instead to 
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the red edge of the adjacent sht. At the assumed po- 
sition, candidate #1 was situated close to a bright star 
that was subtracted before calculating fy. 



5. DISCUSSION 

Are any of our detected objects fluorescent? We know 
from §2 that a good fluorescent candidate should have: i) 
very high EW, ii) possibly (but not necessarily) a double- 
peaked profile, iii) a precise value of its SB. We would also 
like to sec the number density of objects to roughly match 
the density calculated in §2, allowing for the uncertainty 
in clustering. 

As a quick reference, we associated to each candidate 
in Table 3 a "4-" (or "-") if the corresponding indicator 
clearly suggests (or disfavors) a fluorescent origin for the 
Lya emission. Regarding the first criterion we notice 
that, given the sensitivity limits of our V-band image, 
only the two brightest emitters (#1 and #2) have high 
upper limits to their EW (see Table 3) - and it would 
be extremely difficiilt to establish the required EW ob- 
servationally. The low values of the EW for #3 and #6 
indicate a likely non-fluorescent origin for the emission. 

Two objects (#4 and #13) display a very clear signa- 
ture of double-peaked emission. Candidate #4 is com- 
patible with fluorescence if (Tth ^ 60 km s^^. Instead, ob- 
ject #13 is compatible both with fluorescent Lya (with 
(7th ~ 20 km s^^) and Oil emission at z = 0.36. However, 
the high EW, the intensity ratio of the peaks (i.e. that 
the blue peak is brighter than the red peak) all disfavor 
the association with [Oil] emission. 

The comparison of the expected fluorescent SB with 
the estimated SB of the individual candidates is shown 
in Fig. 4. Unfortunately the systemic redshift of the 
quasar is uncertain - z = 3.110 is given by Osmer et al. 
1994, and z = 3.123 by Lanzetta et al. 1993. The ex- 
pected SB from self-shielded clouds (solid line) is calcu- 
lated from Eq. (1), using the corresponding value of L^i, 
(~ lO^^ergs-^Hz-i) and a (0.61 ± 0.02) for our quasar 
(extrapolated from Kuhn et al. 2001). Only the most 
distant objects (4, 5 and 6) have an estimated SB that 
is clearly too high to be due to fluorescence induced by 
the quasar, unless the quasar was brighter in the past. 
Note that two of these (5 and 6) both have /v > Icr (and 
are represented by open-circle symbols in Fig. 4.1), and 
6 in particular has already an EW that disfavors a flu- 
orescent origin. In general, we can notice that the case 
ZQso = 3.110 allows a better fit of the expected relation 
between SB and distance for the detected objects. In 
this case, at least 7—8 candidates have a SB in good 
agreement with the theoretical expectations. 

The theoretical curves shown in Fig. 4.1 are calculated 
assuming isotropic and time-invariant emission from the 
quasar. The candidates beyond the quasar could have 
been illuminated by a different UV flux with respect to 
the value measured by us. In particular, gas clouds at 
20 comoving Mpc from the quasar (within which most 
of our candidates were found) should be responding to 
the ionizing flux of the quasar that was emitted or order 
~ 3 X 10^ yrs before the epoch of observation. At 40 
comoving Mpc, where candidates #4, #5 and #6 lie, the 
time delay is two times larger, i.e. ~ 6 x 10^ yrs, close 
to the expected life-time of the quasar (e.g. Porciani, 



Maghocchetti & Norberg 2004). 

Turning the argument around, the effect of the time 
delay may represent a direct way to measure the radiative 
history of quasars, if the fluorescent origin of the sources 
can be proved. If we assume that our object #4 is truly 
fluorescent, which is plausible given the high EW and the 
clear double-peaked profile, its SB would provide us with 
direct proof that the quasar was shining ^ 10*^ yr ago 
with an ultraviolet luminosity 6-10 times higher than the 
present measurements. Confirmation of the fluorescent 
origin of object #4 would for this reason be especially 
interesting. 

In conclusion, what we can say about the fluorescent 
nature of our detected objects? Unfortunately, there is 
no deflnitive evidence of fluorescence. The EW can be 
poorly constrained for most of our candidates, requiring 
extremely deep broad-band images (for ground observa- 
tions). The SB-distance relation (as well as the emission- 
line profile), given all the theoretical and observational 
uncertainties discussed above, is only an approximate in- 
dicator. 

Nevertheless, collecting all the indication so far, we 
can highlight 5 of the 13 objects which are most prob- 
ably not fluorescent, i.e. #3 (low-EW), #5 (too high 
SB), #6 (both low-EW and too high SB), #11 and #12 
(both objects are unresolved, in front of the quasar and 
with SB at the limit for fiuorescence) . However, given 
their higher upper limit on the EW (mainly because of 
their high line-flux) and their SB, objects #1 and #2 
(if zqso = 3.11) could be good candidates for fluores- 
cence. As noted above, object #4 could also be a good 
candidate, because the profile is flat in the center and 
the presence of a double-peaked profile, compatible with 
fiuorescence if crth ~ 60 km s~^. Finally, objects #7, 
#8, #10, #12 and #13 have SBs compatible with fluo- 
rescent emission (for both quasar redshifts), but the EW 
is poorly constrained. 

As a final remark, it is worth stressing that self- 
shielded fluorescent emitters around a quasar could also 
be detected in absorption as the equivalent of DLAs in 
the field (see §1 for further details). Therefore, the exis- 
tence of such a fluorescent population could in principle 
support the association of DLAs with proto-galactic gas 
clouds (e.g. Haehnelt, Steinmetz & Ranch 1998). How- 
ever, even if all fluorescent emitters around a quasar are 
DLAs, this does not imply the converse, allowing the 
possibility of a different origin for a part of these ab- 
sorbing systems, like disk galaxies already in place (e.g. 
Prochaska & Wolfe 1998). 

5.1. The number density of emitters 

Given our sensitivity limits (SBiim ~ 
10~^^erg s~^cm~^arcsec~^ for a typical aperture 
of 2" X 2") and the assumed value of Lll and a 
for our quasar, we can use the procedure discussed 
in §2 to derive the expected number of fluorescent 
and not-fluorescent sources in our volume. We find 
a physical number density of fiuoresccnt emitters of 
dN/dV ~ 0.3 ± 0.1 Mpc-'\' detectable within Tmax ^ 4.7 
physical Mpc from the quasar. Given the projected 
area covered by our slits in all three configurations 
1.4 physical Mpc^), we therefore should be able 
to detect fluorescent emission within a total volume 
V 13 physical Mpc^. The shape of this volume 
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TABLE 3 

SpECTRO-PHOTOMETRICAL properties of the LYCf CANDIDATES 



Object 




/V±l<7 


EWo 


SB 


Fluorescent?^ 


# 


10~^^ erg s~^cm~'^ 


10~2° erg s-^cm-^A ^ 


A 


10""*^^ erg s~"'^cm~^arcsec~'^ 


EW EP SB 


1 


58.27 ± 6.67 


6.55 ± 5.90 


> 370 


15.81 ± 2.73 


+ + 


2 


18.81 ± 1.20 


3.43 ± 3.69 


> 150 


> 9.72 


+ + ^ 


3 


12.98 ± 1.35 


16.70 ± 3.32 


16t« 


5.80 ± 1.06 


+ 


4 


11.33 ± 1.41 


1.46 ± 4.43 


> 65 


2.68 ± 0.78 


+ + - 


5 


7.30 ± 0.81 


4.76 ± 4.06 


> 36 


2.65 ± 0.49 




6 


6.36 ± 0.68 


8.86 ± 4.06 




1.93 ± 0.40 




7 


5.31 ± 0.82 


0.86 ± 4.06 


>30 


2.71 ± 0.50 


-1- 


8 


4.70 ± 0.67 


0.59 ± 4.43 


>24 


1.84 ± 0.49 


-1- 


9 


4.21 ± 1.16 


0.65 ± 5.17 


>18 


> 2.40 




10 


4.06 ± 1.01 


-5.65 ± 3.69 


>25 


> 2.34 


+ 


11 


3.49 ± 1.05 


-1.98 ± 3.69 


>21 


> 3.52 




12 


2.65 ± 0.56 


3.17 ± 3.32 


>18 


> 1.79 


+ 


13 


2.62 ± 0.45 


2.34 ± 2.58 


>23 {>70<^) 


1.46 ± 0.39 


+ + 


Previous fluorescent hya candidates in literature 


AAl^ 


21.0±5.0 




'^-20 


110.00 ± 26.00 


+ 


FD1<= 


^15.0 




>19 







^in view of the indicators discussed in §2 (wliere EP stays for emission-line profile) 

''if ZQso = 3.110 

'^if the line is Oil at z = 0.36 

'^Adclbcrger ct al. 2006, candidate at 2 = 2.842 
'^Francis & McDonnell 2006, candidate at z = 4.279 



II 


-1 


III 








HI 




Isl 


k 






III 




III 






H ■ 




111 




III 






IID 






III 










I'M 









Fig. 1. — 2D-spectra (in S/N units) of the 13 Lyo candidates (right panels) and the corresponding V-band image of the slit with the 
object positions (white boxes in left panels). The dispersion direction in the spectra is parallel to the x-axis. The scale of the panels is 
32A X 10" for the spectra and 4" X 10" for the images. The white boxes represent the apertures used to calculate the line and V-band 
flux. The values of flux and la noise for each aperture is shown in Table 3 (see text for details). 
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5Q10 



Fig. 2. — V image (left panel), 2D-spectrum (in S/N units, cen- 
tral panel) and calibrated ID-spectrum (solid histogram in right 
panel, dashed lines represent Icr errors) of a bright Oil emitter 
(mv ~ 17) serendipitously found in one of our slits. The emission 
wavelength, the distance between peaks (~ 4A) and the flux ratio 
of the lines, all correspond to Oil doublet at 2 ~ 0.36. The bright 
optical counterpart and the presence of a visible continuum red- 
wards of the double emission line, also indicate that the emitter is 
a low-redshift galaxy. 
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Fig. 3. — ID calibrated spectra of the 13 Lyo candidates. The 

flux (in units of 10~^®erg s^^cm^'^A ) is obtained integrating 
within the aperture represented by the white boxes in the left pan- 
els of Fig.l. The dashed lines represent the Icr noise level. 



is approximatively a cuboid, being constrained from 
the limited field of view of the detector (6'. 8 x 6'. 8, 
corresponding to a physical Mpc scale of 3.2 x 3.2 at this 
redshift). Therefore the number of expected fluorescent 
emitters within r^ax from the quasar is ~ 4 ± 1. Notice, 
however, that this number does not take into account: 
(i) the expected enhancement of the number density of 
objects around the quasar due to clustering effects, (ii) 
any reduction due to beaming of the quasar emission - 
although we note that we would expect the illuminated 
cone to include the line of sight and thus to be more or 
less aligned with the accessible volume. In this context 
we looked to see whether foreground and background 
candidates were located on different sides of the quasar, 
possibly indicating beaming, but found no evidence of 
this. 

By comparison, we expect ^ 6 non-fluorescent sources 
in our total sampled volume if it is similar to the radio 
galaxy environment studied by V05. 

Despite the small statistics and with all the many un- 
certainties, these numbers agree very well with the de- 
tected number of objects in our survey. It may be a 
coincidence, but it is certainly consistent with the indi- 



cation that about half of our detected sources may have 
a fluorescent origin. 



6. SUMMARY 

Based on our recent model of fluorescent Lya emis- 
sion (Paper I), we carried out a Lya survey around the 

2 '--^ 3.1 quasar QSO 0420-388 to search for signatures of 
fluorescence. We used a "multi-slit plus filter" technique 
to sparsely sample a volume of ~ 14000 comoving Mpc"^ 
around the quasar (directly sampling ~ 1700 comoving 
Mpc^). 

We found 13 emission line sources, selected on equiva- 
lent width, which are likely to be Lya at a redshift close 
to the quasar. In order to try to distinguish fluorescent 
objects from internally ionized clouds, we measured three 
possible signatures of fluorescence: i) the line equivalent 
width, ii) the line profile, iii) the SB. We also calculated 
the expected number of fluorescent and non-fluorescent 
sources from theoretical models and recent Lya surveys. 

From a theoretical point of view, the best constraints 
would be a very high EW and, within some limitations, 
the relation of the SB with the distance from the quasar. 
Instead, the line profile cannot give precise information 
on the nature of the sources - a double-peaked profile is a 
possible indication of fiuorescence, but it can be absent, 
or be present in non-fluorescent sources. 

Unfortunately, given the current observational limits, 
there is at present no definitive evidence for fiuores- 
cence. There are good reasons to suspect that 5 ob- 
jects are probably not-fluorescent, and this estimate is in 
good agreement with the number density of such sources 
around radio galaxies. Of the remaining 8 objects, two 
or three of them are quite good fluorescent candidates, 
based on moderately high upper limits to their EW, and, 
in some cases, the presence of extended emission with 
a flat SB and a double-peaked profile. The remaining 
5 — 6 candidates are consistent with fluorescence, but 
their upper-limits on the EW are low enough to provide 
little real constraint. 

On the basis of our model (Paper I) , we were expecting 

3 — 5 fluorescent sources in the inner region around the 
quasar, and, from other surveys of Lya emitters around 
AGNs, we were expecting ~ 6 non-fluorescent sources in 
the total sampled volume. These estimates are certainly 
consistent with the idea that about half of our detected 
sources could be fluorescent. 

In this case, the distance of most of the candidates from 
the ionizing source would imply that the quasar has been 
active for at least 30 Myr, with a UV luminosity similar 
(within a factor of a few) to the present day value. In 
particular, one of the best candidates for fluorescence is 
sufficiently far behind the quasar to imply a life-time of 
> 60 Myr. 

Although there are uncertainties in interpretation, due 
to both observational limitations and theoretical uncer- 
tainties (and no doubt a range of phenomena in Na- 
ture), the current study provides one of the first sta- 
tistical samples of possible fluorescent emitters around 
a high-redshift quasar. Future studies on this category 
of objects will give information about the properties of 
proto-galactic clouds and, furthermore, on the emission 
properties (i.e. spatial and temporal variations) of high- 
redshift quasars. 
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dist from QSO (Comoving Mpc) dist from QSO (Comoving Mpc) 

Fig. 4. — Observed surface brightness (SB) of the 13 Lya candidates against their three-dimensional comoving distance from the quasar 
for two different quasar rcdshifts. The expected SB from self-shielded clouds, as calculated by the theoretical models of Paper I, is shown 
by the solid line. The total errors in the predicted SB should take into account the scatter around the model relation and the uncertainties 
on the derived value of the quasar luminosity (Lll) and spectral shape (a). To give an idea of the effect of these uncertainties in the 
predicted SB we overplot the expected fluorescent SB for a Lyman Limit flux of 2 X Lll (upper dashed line) and 0.5 X Lll (lower dashed 
line). The limiting case of fluorescent emission from a static slab is represented by the dotted line. Open circles indicate the candidates 
with a V-band counterpart with flux greater than Icr for the correspondent photometric aperture (see text for details). The horizontal 
error bars take into account a peculiar velocity shift of ±250km s^^. All the candidates with a SB equal or lower than the theoretical 
expectations curves are compatible with fluorescence induced by the quasar. 
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